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FOREWORD 


This represents the final report on NASA Grant NGL 23-005-336. 

For the most part, the report presents the progress made in the calendar 
year of 1973, although some earlier work is included where needed for 
completeness or continuity. No attempt is made to repeat research findings 
that are already adequately documented in the literature. 

The study was under the direction of Professor J. A. Nicholls, Depart- 
ment of Aerospace Engineering. Dr. R. J. Priem, NASA Lewis Research 
Center, was program manager. 
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ABSTRACT 


This report describes the research conducted in the past year of 
this grant. An experimental study of the passage of a shock wave 
over a burning fuel drop is described. This includes high speed 
framing photographs of the interaction taken at 500, 000 frames per 
second. A theoretical prediction of the ignition of a fuel drop by a 
shock wave is presented and the results compared with earlier 
experimental work. Experimental attempts to generate a detona- 
tion in a liquid fuel drop (kerosene) — liquid oxidizer drop (hydrogen 
peroxide) — inert gas- — environment are described. In this connec- 
tion an appendix is included which gives the analytical prediction of 
power requirements for the drop generator to produce certain size 
drops at a certain mass rate. A Bibliography is also included which 
lists all of the publications resulting from this research grant. 
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I. INTRODUCTION 
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This final report comes at the end of a number of years research 
on the subject of heterogeneous detonations. The motivation for this 
work stemmed from liquid propellant rocket motor combustion instability, 
although certainly the results are also applicable to other jet propulsion 
engines, internal combustion engines, safety aspects of spilled liquid fuel, 
coal mine explosions, airplane crashes, fuel air explosions, etc. The 
research program had been divided into phases and sub-phases, many 
of which had terminated earlier. The topics considered here include the 
passage of a shock wave over a drop which is already burning, a theoret- 
ical treatment of the ignition of a drop by a shock wave, and the experi- 
mental attempt to generate a detonation in a liquid fuel drop — liquid 
oxidizer drop — inert gas mixture. Also an appendix is included which 
presents the theoretical prediction of the power requirements for 
generating liquid drops of a given size and amount. A Bibliography 
is also included which lists all of the publications resulting from this 
research grant. 
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II. RESEARCH RESULTS 


PHASE A. PASSAGE OF A SHOCK WAVE OVER A BURNING DROP 

In rocket motors and other combustion devices, fuel drops are often 
burning when disturbed by a pressure and velocity perturbation. In such 
cases the process of combustion of the drop may be terminated or altered 
as to rate of combustion. Of course, this alters the combustion stability 
characteristics of the motor. In an effort to study this problem an experi- 
mental technique has been developed which allows a single burning fuel 
drop to fall free in the horizontal test section of a shock tube. A controlled 
pressure and velocity pulse is then passed over the drop in the form of a 
shock wave. The ensuing interaction process is then examined through 
use of image converter camera photography, luminosity measurements, 
and (in the past year) ultra high speed schlieren framing camera photography. 
The particular types of interaction studied include the cases of a non- 
ignited drop, the pre -ignited drop which is extinguished by the shock, the 
pre -ignited drop which continues to burn with a wake flame, and the igni- 
tion of a non-burning drop by the shock wave. The results obtained served to 
show the qualitative nature of the interaction, the deformation of the drops, 
the acceleration of the drops, the break-up time, and the threshold condi- 
tions between blow-off and the wake flame conditions. 

A description of the experimental technique and the earlier results 

12 3 

obtained are give n elsewhere 9 3 . However, some of this will be repeated 
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here for the sake of completeness and continuity. The high speed 
photographic study is presented for the first time and includes some 
improvements in the experimental technique. 

1. Experimental Apparatus and Procedure 

A horizontal shock tube consisting of a driver section, driven section, 
dump tank, and diaphragm section was used for these studies. The driver 
section was 3 in. I. D. and 5 ft long which allowed for a sufficient duration 
of uniform convective flow behind the incident shock. The driven section 
was 12 ft long with a rectangular cross section of 1. 5 in. x 2. 5 in. The 
test section was located 9 ft from the end of the diaphragm section with 
schlieren quality plate glass windows 11 3/4 in. long and 1 in. wide. The 
shock tube was fired by venting the gas between the two diaphragms of the 
diaphragm section. 

The procedure for obtaining a time history of deformation and disinte- 
gration of a drop was either taking a series of spark shadow photographs at 
different time intervals after the shock wave intercepted the drop or taking 
a streak photograph through a 1/8 in, wide horizontal slit. Every spark 
photograph taken by the first method was that of a different drop. The run 
to, run reproducibility of the drops was not very good due to the drop injec- 
tion system used but the essential features of the time history of the inter- 
action were retained. Spark photographs were taken with a Beckman and 
Whitley Model 501A image converter camera at the exposure time of 
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1. 0 microsecond. Appropriate timing was obtained by the built-in time 
delay unit of a Tektronix Type 565 dual beam oscilloscope which was 
triggered by the incident shock wave passing over a pressure switch flush 
mounted on the wall of the shock tube upstream. A modified Beckman 
and Whitley Dynafax rotating drum camera was employed for obtaining 
streak pictures. 

The shock speed was determined by the measured time interval for 
the shock wave to travel 2 ft between a pressure switch and a Kistler 
pressure transducer. A pressure switch was used for triggering the 
counting because of its insensitivity to the noise produced by the electric 
spark for the fuel drop ignition. 

For the test of flame blow-off, an RCA 929 photom.ultiplier tube and 
a Texas Instrument LS 400 phototransistor were used. 

The method for producing a suspended ignited drop in the test section 
was based on the sudden retraction of a hypodermic needle from which 
the drop was suspended. The electrode, which produced the ignition spark 
shortly before retraction, was assembled to the needle system and retracted 
with the needle. The needle-electrode assembly is shown in Fig. Land 
the needle tip configuration in Fig. 2. A synchronizing signal from the 
firing of the shock tube activated the ignition spark as well as the solenoid 
which withdrew the locking pin and allowed the needle assembly to retract 
rapidly. Synchronization was accomplished either by the shock wave 
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passing over the pressure switch 8. 05 ft upstream of the drop injection 
position or by the breakage of a thin copper wire on the diaphragm of the 
high pressure section. The method used depended on the test shock 
Mach number. The drop injection system was covered by a dome to 
facilitate evacuation and filling of the shock tube. 

The more recent high speed framing photography study involved some 
changes from the foregoing. In the earlier studies the single drop pro- 
duced often had a tail (thus resembling a tadpole) or the spark discharge 
tended to distort or shatter the drop. Consequently a number of changes 
were made to the drop suspension capillary and electrode. For one thing, 
the piano wire (shown in Fig. 2) was extended beyond the exit of the sus- 
pension capillary. This allowed the drop to be suspended by the wire with 
a consequent marked reduction in wetted area and, hence, deformation to 
the drop when retraction occurred. Also the negative electrode was 
extended to the surface of the drop, thus reducing deformation due to 
spark discharge. A thin mylar shield, painted black, was introduced to 
shield light emitted by the spark from the transient light detector. The 
insulating material used for the positive electrode was replaced by 
teflon to prevent high voltage breakdown. The solid wire used to connect 
the drop suspension capiLary and the retraction springs was replaced by a 
stranded steel wire. This improvement virtually eliminated frequent 
failures of the wire. Finally, the fuel was introduced to the capillary 
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by gravity feed rather than by pressurizing the fuel tank, The net 
result of all of these changes was to produce a good single spherical 
drop, reproducible in size from run to run. 

Experiments involving the use of higher Mach number shocks 
encountered difficulty in synchronizing the firing of the shock tube 
with injection of the fuel drop. If was found that at the higher Mach num- 
bers the shock arrived in the test section before the needle was com- 
pletely retracted, thus resulting in damage to the needle and electrodes. 
Consequently a new technique was developed which uses an exploding 
wire sandwiched between two mylar diaphragms. Energy stored in a 
capacitor (100 joules) is used to explode the wire with consequent rapid 
diaphragm rupture. The sequence of events, then, is as follows. The 
drop was suspended and ignited by a spark at a suitable point in the test 
section and was allowed to fall freely from the capillary. The light 
emitted by the falling -burning -dr op was detected by a transient light 
detector which, in turn, fired the capillary retraction unit and the retraction 
monitor circuit (see Fig. 3). The retraction monitor signal was used to 
activate a delay unit to generate the high energy pulse to rupture the mylar 
diaphragms. A signal from a pressure switch (upstream of the test section) 
was used to trigger the oscilloscope and delayed output from the scope was 
used to trigger a pulse light output from the xenon lamp. 
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For the high speed photographic study a Beckman-Whitley 330 framing 
camera was used at a speed of 480, 000 frames/ sec for an approximate 
duration of 160 jisec in order to record the interaction of the shock with 
the drop. Schlieren optics were employed. Of course, this allowed the 
complete interaction to be obtained in one run, whereas with the earlier 
image camera records only one photograph per run could be made. 

For the experiments, the test section gas was either oxygen, nitrogen, 

or air, with various initial pressure levels, P^, possible. Various shock 

Mach numbers, M , were generated. The fuel used was diethylcyclo- 

s 

hexane (DECK). 

2. Results and Discussion 

Image converter camera photographs of the interaction of an 
M = 2. 02 shock with a 1. 6 mm DECK drop in one atmosphere of 0 0 
are shown in Fig. 4. A weak bow shock can be seen in some of the 
photographs. Disturbances from secondary drops, or the "tail M of the 
parent drop when the needle is suddenly retracted, are to be noted in 
some cases. It is noted that the incident shock wave is accelerated in 
the hot region (see Run Nos. 31, 20 and 55). The initial flame appears 
to be laminar but becomes turbulent as it is blown off the parent drop and 
into the wake by the convective flow behind the shock wave. The initial 
flame apparently ignites the cloud of smaller drops shed from the parent 
drop (micromist)(see Run No. 55). The flame attached to the wake of the 
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disintegrating drop (wake flame) becomes fully developed when the 
rate of micromist shedding attains some level, as is seen in Run Nos. 

53 and 25. The parent drop apparently acts as a flame holder. At 
the later stages of interaction, the bow shock is reflected from the 
bottom and the top of the shock tube walls as is seen in the pictures 
of Run Nos. 24, 25, and 53. Comparison of these results with similar 
ones for non-ignited drops shows that the lateral deformation of drops 
is larger in initially ignited drops than non-ignited ones, and that the 
micromist shed from the drops which have a wake flame is consumed 
much faster than that from non-ignited drops. 

Under certain conditions of drop size, shock Mach number, and 
initial pressure, the flame will not stabilize in the wake but will be 
blown off. In those cases the disintegrating drop rapidly takes on the 
appearance of the non-ignited drop. 

Streak photographs show a continuous recording of the interaction 
of a shock with one drop. Figures 5 and 6 are streak schlieren photo- 
graphs of non-ignited and ignited drops in pure oxygen with an initial 
pressure of 30 in. Hg. The shock Mach number is about 1. 7 so that 
the flow behind the shock wave is subsonic and hence compression waves 
are seen to be transmitted upstream. Comparison of Figs. 5 and 6 
shows an apparent effect of the wake flame on break-up characteristics. 
The curvature of the trajectory of the drop at the leading edge changes 
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abruptly at about the end of the break-up on the streak photograph or on 
the x-t plane. The break-up time of a drop, t^, which is defined as the 
time from the instant when the incident shock intercepts the drop until 
the disintegrating drop is accelerated to 60% of the convective flow veloc- 
ity behind the incident shock wave seems to be longer in the ease of a 
drop which has a wake flame. The opaqueness of the wake in Fig. 5 is 
attributable to the microspray droplets whereas in Fig. 6 these drops 
have vaporized and burned. 

Figure 7 shows the deformation of drops (i. e. , the ratio of the maxi- 
mum width of dense area d perpendicular to the flow to the initial diam- 
eter dg) with time elapsed from shock passage through the center of the 
drop. Data is shown for drops with a wake flame as well as for non~ 

ignited drops. The solid lines in the figure are the correlation of Ranger 
§ 

and Nicholls for water drops. The deformation of the non-ignited drops 
are in good agreement with this correlation but, as noted, the drops with 
a wake flame produce much larger deformations. Similar data obtained 
for the "blow-off 1 case indicated somewhat greater deformation than for 
the non-ignited case but not as large as that of the drops with a wake 
flame. The surface tension of DECK is 28 dynes/ cm, while that of water 
is 73 dynes/cm at 20°C. In view of the fact that the DECH deformation 
is correlated well with the equation derived from experiment with water, 
the effect of surface tension may not be large. Also, the viscosities of 
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DECH and "water are about the same. Hence, it is tentatively concluded 
that the differences in deformation are largely attributable to altered 
pressure distributions around the drop as a consequence of the flame. 

In spite of the fact that the original drop is changing shape and losing 
mass by micromist shedding, it is interesting to consider the drop as a 
spherical particle of constant mass. Then the displacement of the drop 
can be expressed as 



or, x = KT 

where, x = dimensionless distance 
T = dimensionless time 

K = 3/8 , and is the apparent drag coefficient based on the 

D o D o 

initial diameter of the drop, x is the displacement of the leading edge 
of the drop along the flow from the original position, and p 2 are the 
velocity and density behind the incident shock, respectively, d Q is the 
original drop diameter, p^ is the density of the liquid drop, and t is the 
time after the beginning of interaction. Constant acceleration is assumed 
and relative velocity between gas and drop, drag coefficient, and drop 
diameter are evaluated at t - 0. 
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The displacement of the drop satisfies Eq. (1) under all the test 

conditions reported here for a good fraction of the break-up time. 

Figure 8 is an example. The apparent drag coefficient was determined 

from such displacement data and indicated little variation with Reynolds 

number. The Mach number of the convective flow did vary some in 

various runs but this effect was not isolated. The averaged values of 

C n are: 1. 80 for drops with a wake flame, 2. 10 for non-ignited drops, 

0 

and 2. 66 for drops with "blown off" flames. The drag coefficient for 
the drops with a wake flame is the smallest, probably due to the altered 
pressure distribution. 

Drop break-up time, t^, was found to be somewhat larger for drops 
with a wake flame. 

For the determination of "blow off" limit, two types of photoelectric 
sensors were used. An RCA 929 photomultiplier tube viewed the whole 
window of the test section and a Texas Instrument LS 400 phototransistor 
viewed a small portion at the original drop position. It was easy to 
conclude from the voltage time outputs of these sensors whether a wake 
flame or blow-off was realized. 

Figure 9 shows the line of demarcation between the wake flame and 
blow-off regions in the initial pressure — shock Mach number plane. The 
limiting initial pressure increases at about ~ 2, where the convective 
flow behind the incident shock is sonic. Otherwise, the limiting initial 
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pressure is fairly insensitive to the Mach number. Other experiments, 
conducted with a 75% 0,,-25% mixture, yielded very similar behavior 
but with higher limiting initial pressure values. 

The high speed framing camera results will now be discussed. The 
various cases studied are indicated in Table I. 

Table I. Test Conditions, High Speed Photography 


Case 

Initial 

Condition 

Test 

Section 

Medium 

Mach 

No. 

Pressure 
In. /Hg 

After 

Interaction 

Drop Size 
/im 

(approx. ) 

a. 

Non-ignited 

Nitrogen 

2.05 

29.1 

Non-ignited 

1500 

b. 

Ignited 

Air 

2.05 

29.1 

Extinguished 1250 

c. 

Ignited 

Oxygen 

2.05 

29.1 

Burning 

1250 

d. 

Non-ignited 

Nitrogen 

3. 5 

11. 5 

Non-ignited 

1500 

e. 

Non-ignited 

Oxygen 

3. 5 

11. 5 

Ignited 

1500 


Photographic sequences of these various cases were obtained. 

Figure 10 shows the comparison between cases a, b, and c, each of 
them taken at the same time after interacting with the shock ( ! 0 T micro- 
seconds corresponds to the instant when the shock hits the parent drop). 

In these photographs there is, for the most part, only one drop. However, 
the window quality is not good and there are smudges and vertical marker 
lines on the windows. The pre-ignited drop is injected approximately 
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1/ 2 in. above the center line of the test section and burns for about 
100 milliseconds before interaction with the shock. This is the reason 
for the presence of column of hot gases in the cases b and c. When the 
incident shock enters the region of hot gases surrounding the drop, 
depending on the acoustic impedance of the cold uniform region and the 
hot gases in the flame, transmitted and reflected waves will be obtained 
at the interface. Hence, it is plausible to expect a significant difference 
between the wake structures in the cases a, b, and c, as can be seen in 
Fig. 10. The upward shift in the wake of the drop is not because of buoy- 
ancy effect (since Grashoff number « 1) but due to the altered shock 
structure and shape when it enters the hot region. The combustion 
products generated before shock arrival are seen to blow downstream. 

In case b the flame is blown-off and no further combustion occurs; the 
subsequent drop breakup process then resembles case a. In case c 
combustion continues in the wake of the drop. 

The structure of the wake is considerably affected by the drop break- 
up characteristics and also by the extent of micromist evaporation. In 
the case of the wake burning in oxygen, the micromist evaporation 
process will be enhanced due to higher flame temperature and thus micro- 
mist will contribute to the mass of gaseous contents of the recirculation, 
thereby reducing the mass entrained from the external convective flow. 
This seems to be the reason for a relatively longer and wider wake 
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structure in case c. The resultant shallower wake angle can be expected 
to affect the pressure distribution in the wake considerably because of 
its effect on the deflection of the external convective flow. In turn, 
pressure distribution can affect drop acceleration and drop breakup 
characteristics. 

Cases d and e are shown in Fig. 11 wherein neither drop was burning 
initially. The strong shock was sufficient to ignite the drop when in O^, 
case e. Ignition was determined by an appreciable photocell output. In 
the photographs, ignition is not obvious. However, it must have occurred 
somewhat before 80 (xs in that the length of the visible wake is much less 
in the last photograph and the shape of the bow shock has been changed. 
There is some indication of a blast wave interacting with the bow shock 
in this photograph. Photographs taken at later times indicate the much 
faster consumption of the drop in the burning case than in the case. 

The differences in the breakup pattern for non-burning drops for the 
two Mach numbers can be noted by comparing Figs. 10 and 11 for the cases 
a and d. In both cases it is possible to see the formation of surface waves. 
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PHASE B. ENERGY RELEASE PATTERNS 


1. Theory of Shock Ignition, of Fuel Drops 
A. Introduction 

The essential nature of two-phase detonations is ultimately connected 

with the manner in which the liquid phase reactant is converted to its vapor 

phase, and to how the two vapor phase reactants then mix and ignite. In 

the case of two-phase detonations in which the liquid phase reactant is 

initially in the form of a spray (of fuel droplets, say), the conversion process 

involved is that of the interaction of the shock wave with each of the drops 

individually in the spray. Much empirical data have been gathered which 

have shown the major characteristics of this particular process for the case 

4-8 

of shock interaction with nom-reacting drops . Unfortunately, must less 
information is available for the case of reacting drops, which forms the 
subject of interest her e^ 

Both non-reacting and reacting drops share certain features of a shock- 
wave interaction. As the shock wave is passing over the initially stationary 
droplet, very little of consequence occurs. In pai'tieular, the drop does not 
change its velocity during this period while the surrounding gas is acceler- 
ated. This produces a flow field around the droplet. If the shock strength 
is sufficient, the free stream velocity can be supersonic with respect to 
the drop. Under these circumstances, a bow shock, wake shock, and other 
characteristics of supersonic flow over a sphere are apparent (Fig. 12). 
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The drop subsequently contracts along an axis which is parallel to the free 
stream, and expands in the transverse direction, so that the frontal area 
exposed to the flow increases with time. The drop also begins to accelerate 
in the direction of the moving shock. 

The fluid near the surface in the drop is set in motion by the boundary 
layer gases of the convective flow. The liquid boundary layer separates 
from the droplet and is apparently broken up into a spray of droplets whose 
diameter is of the order of the boundary layer thickness. This "microspray" 
{(i spray) is thereby introduced into the near wake region of the parent drop- 
let. It appears that, at the same time, a Taylor instability occurs on the 

front surface of the parent drop, producing waves whose amplitudes grow 

13 

and will eventually cause "catastrophic" disintegration of the parent drop . 

After a certain period of time following passage of the incident shock, 
the parent drop ceases its transverse growth; the frontal diameter there- 
after decreases. Whether or not catastrophic breakup due to Taylor insta- 
bilities has occurred by this point is not clear. Nonetheless, if the parent 
drop is reactive, it is somewhat after this time that the ignition of evaporated 
ix spray occurs in the near wake region. The ignition is explosive in charac- 
ter, producing a blast wave. 

Beyond this ignition, the process of mass stripping from the parent 
drop (or its fragments in the event of catastrophic breakup) continues until 
nothing remains of it. If ignition has occurred, in the reactive case, the 
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continued stripping supplies fresh fuel to the wake region where it is then 
consumed. This post-ignition combustion process is usually relatively 
"smooth’*; occasionally, however, multiple explosive ignitions occur in 
sequence. 

Only incident shocks of sufficient strength produce explosive ignition. 

As a rule, M > 3 is required, typically. In what follows, it will there- 
s 

fore be assumed that the incident shock is of this strength, and also that 
the resulting flow field around the parent droplet is initially supersonic. 

The capacity to describe the sequence of events, which was qualita- 
tively outlined above, in sufficient detail as to allow reasonable prediction 

of the time between incident shock contact with a reactive drop and its 

14 

explosive wake ignition, is of obvious interest. Kauffman suggested a 
model in which groups of stripped juspray, which move rearward into the 
wake region, remain as they move in volume elements of fixed identity and 
fixed geometry. The n spray in each such volume element evaporates, 
and the vapor produced mixes homogeneously throughout the element. 
Ignition is identified with that element in which the concentration of fuel 
is largest; ignition times are assessed from the interval separating inci- 
dent shock contact and the formation of this element. 

15 

Pierce treated a simple model in which stripped p spray enters the 
wake and therein evaporates and reacts with the hot oxidizer. The energy 
thus liberated is conducted away to the external flow until the total evapora 
tion rate (which controls the reaction rate) becomes so large that thermal 
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energy begins to accumulate in the wake region. The point in time at which 
this occurs is described to be the instant of ignition, since as the tempera- 
ture of the region increases, the processes responsible for further energy 
liberation are accelerated. 

Recently, Fishburn has suggested that the density of jn spray drops 
(in the wake region of the parent drop) during the first period following inci- 
dent shock passage, produces local fuel vapor concentrations which are 
too low to support explosive ignition. That is, boundary layer stripping of 
the contiguous parent droplet is considered to be too slow a process to 
allow for high y . spray densities. Therefore it is argued that fragmenta- 
tion of the parent drop, through Taylor instabilities, mast first occur; each 
fragment then individually strips, resulting in a much higher rate of ju spray 
production. Again, the ft spray drops enter the wake, evaporate and react, 
and a chemical induction time is added (to the time until fragmentation 
occurs) to obtain the overall ignition delay. 

Each of these theories shows some degree of agreement with the avail - 

9 11 14 

able experimental ignition time data } 9 , yet none, in itself, is able to 

convincingly explain all of the observations made of the ignition process. 
The ability to extend any one of these computations to reliably include a 
wide range of fuel/oxidizer combinations is rather questionable. It is the 
purpose of the present exploratory study to contribute towards the develop- 
ment of a unified theory of shock induced reactive droplet ignition. 


IS 


B. Theoretical Model 


In its present form the sequence of events leading to explosive ignition 
is considered as comprised of two separate intervals; n, tely, a "dormant" 
period, followed by an "active" period. As such, this model is similar in 
kind to Fishburn's two-stage model. That is, if the same terminology were 
applied to that model, the "dormant" and "active" periods would correspond 
to the intervals before and after parent drop fragmentation, respectively. 

It is to be emphasized at the outset that the two periods in the following 
formulation both differ qualitatively (as well as quantitatively) from the 
Fishburn model. 

(i) Dormant Period 

This period begins upon initial contact by the incident shock with the 
spherical parent droplet in its undisturbed, motionless state. Four pro- 
cesses are initiated. First, the droplet begins to accelerate. This in 
turn initiates the development of Taylor instabilities on the forward sur- 
face. At the same time, the droplet begins to flatten, and the boundary 
layer formation in the liquid surface commences. 

Complete boundary layer formation requires a nontrivial induction 
17 

time , however, it is believed that the mass stripping process begins well 
before the liquid boundary layer is fully developed. This is based on the 
many early-time photographs such as in Fig. 13, as well as on mass loss 

Q 

measurements . In any event, it is clear that the mass loss rate acceler- 

g 

ates with increasing time during this period. In fact, Reinecke has obtained 
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a reasonabx 


r i.: ;rical correlation for mass loss, assessed from x-ray 
photographs of stripping water drops, which is 

ir = I [i + cos . (t/t s )] . (1) 

0 

The mass, m, of the droplet at any time, t, during the breakup process 

is thus correlated with its initial mass, m , and the time to complete 

disintegration, t (stripping time), 
s 

As the droplet continues to flatten, its frontal diameter increases 
rapidly. Fluid in the liquid boundary layer, which travels from the for- 
ward stagnation point to the maximum perimeter before separating into 
ju spray, thus travels progressively further before being stripped off. 

This, as well as the fact that the boundary layer becomes more fully 
developed with passing time, leads to the expectation of increased pspray 
size with increasing time. 

No widely accepted means of computing /i spray sizes and separation 

velocities exists. There is some agreement, however, that the juspray 

diameter should be of the order of the liquid boundary layer thickness 

14 , 15 

just prior to its separation * . Most analytical estimates predict /ispray 

diameters of the order of .01 of the parent drop diameter, which roughly 

4 7 

agrees with what experimental evidence is available * . 

For example, the separation process can be envisioned to proceed in 
two stages: First, an annular sheet or film of fluid is shed from the drop 


20 



periphery, and second, this sheet breaks up into the p spray drops by 

means of a process similar to that which occurs in the breakup of a free 
18 

liquid jet . The micromist drops can, by that analogy, be expected to 
have diameters which are approximately twice the thickness of the sheet. 
The sheet thickness itself is estimated by imposing conservation of mass 
and momentum on the fluid which enters it, between the positions just 
before separation (when the fluid is in the parent drop boundary layer) 
and just after separation (when the fluid is in the annular sheet and has 
a flat velocity profile) . 

When the liquid boundary layer is assumed to have a Taylor velocity 
5 

profile , this particular means of analysis produces reasonable early- 

15 

time mass loss rates and jaspray sizes . In addition, the velocity of the 
spray when it separates from the parent drop, u^, can be readily derived 
to the form 


U d A 

u c = 72 


( 2 ) 


In which u is the relative convective flow velocity, while 

w 


1/3 

(W\ 


( 3 ) 


and Pp p„, fig are the densities and viscosities of the liquid and of the 
free stream convective flow. 
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It is reasonably clear that boundary layer stripping cannot account for 

the mass loss rates over the entire breakup time. The surface wave con- 

19 13 

cept of Collins or the Taylor instability theory of Fishburn is needed 

at later times to explain the high stripping rates. It is believed, however, 

that simple boundary layer stripping predominates the early stages of 

breakup, and Eq. (2) , in spite of its crudeness, serves to show that 

u,/u =0(10 for typical intermediate shock strengths and typical fuel/ 

Cl c 

oxidizer combinations. That is, the fluid separates from the parent drop 
at velocities which are always much lower than the prevailing free stream 
gas velocity. 

The velocity at which the point of separation moves in the transverse 
(or, radial) direction due to drop flattening may be evaluated from the well, 
accepted empirical form^ 

D = 1+kT (4) 

in which D is the ratio of parent drop frontal diameter at time t to its 
initial (undisturbed) diameter, D q , and 

T=J^ 2 t (5) 

0 


is the non*-dimensional time, J3 = P^/Pj)* is the free stream gas velocity 
relative to the drop at t - 0, and a is a correlating coefficient whose value 
is a ~ 1. 70. From Eq. (4), the radial velocity, u , of the separation point is 
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( 6 ) 


u r /u c = (1/2) otj3 ly/2 (u 2 /u c ) . 

The radial velocities computed from Eq. (6) show that u /u is also 

!T C 

0(10 “V 

The inference from these two simple results is that in general the 
separated liquid film can be expected to turn in the direction of the pre- 
vailing local flow field while it is still contiguous with the parent drop 
(i. e. before breaking up into juspray drops) . When the external flow field 
is subsonic, the /xspray is in this fashion carried rearward, more or less 
parallel to the axis of symmetry. It is gradually accelerated, and ulti- 
mately reaches the convective flow velocity. An example of this is shown 
in Fig. 14. 

When the flow field is supersonic, which is the case of interest here, 
it appears that the separated film is turned inward by the flow structure, 
in such a way that the p;spray is initially carried into the free shear layer 
of the near wake, above the recirculation zone. It subsequently becomes 
engulfed in the expanding recirculation zone itself, and most does not 
escape, having entered at low velocity. This results in a remarkably well- 
defined recirculation zone, as in Figs. 15 and 15. 

In fact, the rearward velocity of the recirculation zone tip can be 
estimated by noting that its geometry remains roughly constant while the 
parent drop is growing. Since the droplet growth rate is given by Eq. (6), 
the rate of growth in recirculation zone length is simply 
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dL/dt = u^ ctn j3 


n 


in which J3 is the recirculation zone angle with respect to the axis. For 
typical supersonic wakes, 4 < ctn |3 < 6. Equation (7) then shows that 
dL/dt <K u^, even for the larger value of ctn 0. The actual velocity of 
the recirculation zone tip, measured from streak photographs such as 
Fig. 15(b), agrees with the order of magnitude predicted by Eq. (7) . (It 
should be observed that the recirculation zone is not filled by microspray in 
the subsonic case and so it is not visible on photographs such as Fig. 14. 

Within the recirculation zone, the trapped jr spray move with the low 
velocity vortices and evaporate. If a given pspray drop is formed at time r 
(measured. after incident shock contact with the parent drop), with initial 
diameter dM (t), then its diameter at time t (i. e. after an interval t-r) 


can be approximated by the quiescent evaporation rate expression 

1 1/2 
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i 


d M( 2(r) - 2C (t 


-r)] : 


0 


t < t 

— ex 

t> t 

— ex 


( 8 ) 


in which 


C = -^— In 
C p p £ 


[l + J(T g -T £ )] , 


(9) 


where k is the coefficient of thermal conductivity of the gas, T is its 

& 

temperature, and its specific heat, and £ and are the latent heat 
and temperature of the liquid pspray drop. The time to extinction of the 


pspray drop (complete evaporation), t , is 

ex 
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( 10 ) 


t 

ex 


- T = 


d M 0 2 < T ) 

2C 


The rate of mass evaporation from this p spray drop at time t is 

m ey (t, r) = C 2 d M (t, r) , (11) 


where Cg - irp^C/2. 

These formulae allow computation of the total amount of fuel vapor which 
has evaporated from all jul spray drops in the recirculation zone up to time t, 
taking into account the time varying initial pspray size, in the following man- 
ner. First, the number of pspray drops, fin, which are formed during an 
interval, 6r, about time r, is 


5n(r) = 


m(r) 5 t 

m M W 

o 


( 12 ) 


where m(r) is the parent drop stripping rate, at t = t, and m^ (r) is the 
mass of the p spray drops formed at that time; i. e. mjj (r) - p^ ir d]vi 0 (r)/6. 
The contribution at time t > r to the total pspray evaporation rate from 
p spray drops formed at t = r is 


or, in the limit, 


5m (t) a m (t,r) 5n(r) 

ev ev 


dm (t) « CL d, - (t, r) ■ dr 

ev w 2 M v (r) 

o 


( 13 ) 
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Therefore, the total evaporation rate, at time t, due to all ju spray drops 
present in the recirculation zone at that time, is found by integrating Eq. 
(18) over all r t; namely, 


m (t) = C 9 
ev 2 


L 

/ 

T— 0 


m(r) 
m M (T) 

IYIq 


dT 


( 14 ) 


and the total mass of fuel vapor evaporated from all p spray drops prior to 
time t is 


m (t) = Co 
ev 2 


t t* 

// 

t*=G t=0 


m(r) d M (t*,r) 


dr dt* 


(15) 


A simpl.? parent drop stripping rate which can be used is 


jtm 

m (r) = sm ~ > 
s s 


(16) 


which derives from Eq. (1), by differentiation, and the time varying p. spray 
size, dj^ (r) is obtained in the fashion described earlier. Then, upon 
appropriate selection of values for the parameters which appear in Eq. (9), 
the vapor accumulation in the recirculation zone, m eY , can be computed 
from Eq. (15) as a function of time. 

It is noted parenthetically that mj^ is, of course, never actually zero, 
but that, on the other hand, m(r) is zero at r = 0, That is, actual parent 
drop stripping can begin only when the liquid boundary layer has developed 



to the point that Its kinetic energy cannot be dissipated by the liquid surface 
at the point of separation. This occurs, approximately, when the balance 

d 

s 

So-jj ir D = (1/2) p^ttD f u £ 2 (y, t) dy (17) 

0 

occurs, in which is the liquid surface tension, u^fy, t) is the time-varying 

liquid boundary layer velocity profile (measured inward from the surface 

at the point of separation), D is the parent drop frontal diameter, and d is 

s 

the thickness of the annular sheet, when it is separated. This concept of a 

g 

"boundary layer induction time" is indeed similar in spirit to that of Ranger , 
but it d^es not require that the boundary layer be fully established before 
initiation of stripping. For the present purposes, an estimate of the mini- 
mum film thickness, from Eq. (17), was used to compute Kij^ o (0) required 
in Eq. (15). 

Characteristic calculations from Eq. (16) are summarized in Figs. 17 
and 18, for the case of diethyl cyclohexane (DECH) drops in oxygen. The 
vapor accumulation, at the moment, t*, of maximum droplet expansion, is 
shown in Fig. 17 as a fraction of the total mass removed from the parent 
drop up to that time, Figure 18 shows the maximum p spray initial 

diameter, dj£ Q (t*). The effect of initial pressure, F^, as well as parent 
drop initU diameter, D 0 , and incident shock Mach number, M is 
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demonstrated. The extent of the evaporation which has occurred by t* 
is surprisingly small for D q somewhat larger than 300jut. 

The equivalence ratio, <f>, corresponding to the accumulation of fuel 
vapor within the recirculation zone at t*, is shown on Fig. 19. Equiva- 
lence ratio is defined as the quotient of actual fuel/oxidizer ratio to the 
stoichiometric fuel/oxidizer ratio, so that = 1 represents the stoichio- 
metric condition. It is observed that at reduced pressures and for small 
parent drop sizes, the equivalence ratio could reach significant values. 

The calculations indicate, however, that the accumulation rates are slow; 
i. e. m r (t) from Eq. (14) is never very large for parent drops which are 
larger than 300p diameter. 

From these results, three conclusions can be drawn. First, the 
recirculation zone does not impulsively reach an extremely fuel rich condi- 
tion (as was initially suspected) . On the contrary, the region is generally 
quite lean; second, in those cases for which m is non-negligible, con- 
tinuous reaction with the hot oxidizer will preclude accumulation of fuel 
vapor to the extent indicated in Fig. 19. Moreover, at time t*, most 
of the liquid fuel which was stripped from the parent drop resides in its 
recirculation zone in the form of unevaporated jut spray. Therefore, there is 
little possibility that the recirculation zone can itself support an explosive 
ignition. In fact, when ignition occurs, the entire wake region is consumed 
except for the recirculation zone, as can be seen in Fig. 20. 
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(ii) Active Period 

The frontal diameter of the parent drop reaches its maximum when 

t = t*; t* correlates roughly with T - 1. 3. Catastrophic breakup, as pre- 

13 

dieted by Fishburn , would have occurred prior to this time; viz. , approxi- 
mately at T = 0, 9. It is not completely clear as to whether or not this is 
in fact the case, and, in any event, the aggregate representing the remain- 
ing parent drop shows no abrupt change in velocity during this period. If 
it has been shattered into fragments, these fragments, as a closely packed 
group, evidently, thereafter, behave as a single (porous) body. Fragmen- 
tation would, of course, explain the accelerated stripping rates (which occur 

19 

at about this time), but Collins’ surface wave theory does so equally well. 
The experimental data available at this writing simply do not allow for dis- 
crimination between these two possibilities. 

Nevertheless, it is not essential to the present phenomenological 
description of ignition that the cause for accelerated mass removal be actually 
identified, because a change in the stripping mechanism does not appear to 
be the single event that is primarily responsible for eventual explosive igni- 
tion. Instead, it appears that the termination of parent drop flattening, 
which allows the escape of substantial quantities of p spray from the near 
wake region, is responsible for the final events leading to ignition. Thdt is, 
p spray which is shed during the period following time t* is not engulfed by 
an expanding recirculation zone, but rather is injected outside of a shrinking 
near wake so that it becomes exposed to the high speed convective flow. 
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As the escaping fi spray begins to move rearward, it is itself acceler- 
ated, and this requires a small, but finite, amount of time. With a drag 
coefficient of unity, the equation of motion for a ju spray drop is simply 


du 

c 

dt 


a 

4 P i d M 


(18) 


in terms of the relative convective velocity, u c * For the case of constant 
this integrates to 


u 

c 




u t \ 1 

3 p 2 C o \ 

4p l d M / 


( 19 ) 


where u c is the relative velocity at separation. Equation (19) provides 

G 

an upper bound on acceleration times. For u Cq = u^, the time for a lOp 

DECH drop to reach u /u = 0. 3 is 0(1 jusec) corresponding to M =4, 

c c s 

o 

and = 1 atm oxygen; the time to reach u c / u Cq = 0.1 is 0(5 jusec). 
Acceleration times scale linearly with ju spray diameter. 

The {i spray that escapes at t > t* can be observed on streak photo- 
graphs such as Fig. 15(b). Although the acceleration time of an escaped 
/xspray drop is very short, it is largely converted to the vapor phase dur- 
ing that period. The 10 fi DECH drop in the above example has a Weber 
number. We = (p^ u ^ which is 0(10 ^ and this is still far above 

the minimum Weber number that corresponds to the stripping mode, 
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We ~ 15. If the stripping mechanisms are not essentially different than those 

for larger drops, the time for this lOpdrop to strip is 0(1 psec). This is 

the same time order as its acceleration period. The "second generation" 

pspray produced during this stripping process, whose diameters should be 
-I 

0(10 p), would vaporize in negligiule time orders. 

If the stripping mechanism is unacceptable for these small drop sizes, 

rapid conversion of the escaped pspray to the vapor phase can also be 

explained by convection assisted evaporation during the acceleration period. 

When the relative convective velocity is appreciable, as in the case of 

escaped pspray drops, the quiescent evaporation rate, given by Eq. (8), 

22 

should be replaced with 

^ = C (1+Re l/2 pr l/3, > < 20) 

in which Re = (p^ u c d-^/pg, and /k. The second term in brackets 

can account for a five-fold increase in evaporation rate under typical condi- 
tions, reducing the characteristic evaporation time of the lOp diameter drop 
from 200 psec to 40 psec. Hence, if evaporation alone must account for 
the conversion of the p spray to vapor, an appreciable amount of this conver- 
sion can be shown to occur during the acceleration period. 

It is most probable that both stripping and evaporation occur simultan- 
eously. Neither appears to cause complete conversion of the p spray to vapor; 
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that is, the jn spray remains visible on photographs. However, for modeling 
purposes, complete conversion will be assumed. 

Under that assumption, the essential feature of the parent drop disinte- 
gration process, beginning at t*, is characterized by the impulsive continual 
injection of relatively large quantities of reactive vapor into the outer near 
wake region of the parent drop. This is somewhat similar to the occurrence 

of forward stagnation point mass addition as it might occur in a supersonic 

23 

flow about a solid sphere or cylinder . By that comparison, vapor concen- 
trations can be expected to decrease with increasing axial and radial dis- 
tance within the wake; the greater variation would be in the radial direction. 
In fact, the assumption of constant radial vapor di stribution with increasing 
axial distance is not wholly unreasonable, over the first few diameters of 
length. On the other hand, an assumption of radial uniformity in the vapor 
concentration would appear to be a rather dangerous oversimplification. 

Each local element of mixed fuel vapor and oxidizer may for modeling 
purposes be regarded as a homogeneous chemical system which, at the 
instant of initial mixing, has a specific mixture ratio and initial tempera- 
ture. As reaction proceeds toward ’’equilibrium’’, the products remain 
within the element. Diffusion and thermal conduction between adjacent ele- 
ments are ignored. The later was justified based on an estimate of laminar 
heat transfer over the time orders of interest. Each element therefore 
experiences an accelerated reaction rate due to self -heating. After an 
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induction time, the maximum rate of temperature rise is reached. 

The energy release rate at this time can be of explosive proportions; if so, 
it marks the point of wake ignition. 

Now, the initial temperature in each element is a function of the vapor 
concentration there. That is, the maximum temperature of the fuel in the 
condensed phase is its boiling temperature at the prevailing local pressure. 
When the fuel changes phase and mixes with the oxidizer, the mixture 
temperature (assuming a mixing process which takes place at constant 
pressure) can readily be shown to be 


T T b - ^ C pf , gT 2 

X 1 +0! 1 + 0t 


( 21 ) 


in which T^ is the boiling point of the liquid, T^ Is the static temperature 
of the oxidizer in the free stream prior to mixing, and 



Here, k £ is the mass fraction of fuel vapor in the element, and C , C „ 
f J px pf 

are the constant pressure specific heats of the fuel vapor and oxidizer, 
respectively. 

The self-heating process is initiated at this temperature. To assess 

the duration of the induction period, the Edelman -Fortune global reaction 

24 

rate equation for hydrocarbon vapor combinations with oxygen is used. 
This is 
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9Cp 

i 

"aT 


chm 


c p- <5 “.825 ^ ~l/2 ^ 

5, 52 x 10 p T C f Cq exp 

2 


[(- 1.22x IoVt] » (23) 


in which and are the fuel and oxidizer concentrations in gram -moles/ 
3 ^ 

cm , while the units of system pressure and temperature are atmospheres 
and Kelvin degrees. 

The induction period is assumed to proceed under a constant pressure 
condition. Thus, as heat is liberated by reaction, the element volume in- 
creases. Concentration changes are then due both to volumetric expansion 
and chemical reaction; the rate of energy liberation is determined from the 
latter. The fuel concentration in a constant pressure system can not be 
determined as a function of time from Eq. (23) alone. However, T(t) may 
be obtained from this expression directly, because p - p(T) , where p is the 
total element mass density. 

For any hydrocarbon, the stoichiometric equation is 

C n H m + (n + 1/4 m) C> 2 - n C0 2 + 1/2 m HgO , 

and, temporarily assuming a fuel-lean system, the disappearance rate 
of oxygen molecules is assumed to be approximately governed by this 
equation. Hence, 


at 


= n 


chm 


0C, 

i 

s at 


chm 


(24) 
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in which n g = n -f m/4. Now for either species reactant, C = p./w.., 


where p^ is its mass density and w^ its molecular weight. In terms of 
its mass fraction jc, this is C. = /n p/w^. Hence, Eq. (24) becomes 


3k 


O, 


w. 


at 


V K f 


. w- at 

chm f 


chm 


( 25 ) 


and also Eq. (23) may be rewritten in the form 


1 3K f 


w f at 


chm 


1/2 

K f K a 

R in 8 -,825 1/2 U 2 

= - 5. 52 x 10 p T p 


1/2 

w w 
f O t 


♦ exp [(- 1. 22 x 10^)/T] . 


(26) 


Of course, mass fractions are not affected by simple volumetric expansion, 
and so Eq. (25) can be integrated to give 


K 0 ~ K 0 


W 

„ " ~vT. (K f. ' K i! ■ 

2. f i 

1 


(27) 


Note that in Eq, (26), the units of p are gram/cni . 

Now, assuming a fuel lean condition, the heat of combustion per unit 

mass of fuel, 3/ , is approximately constant. Temperature increases due 
c 

to combustion are then related to changes in quantity of fuel present through 
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d/c, 


= '"d 2, 

chm /7 c 


dT 


(28) 


which is integrated to give 


* £ -*l -# (T “ T ^ > 


(29) 


i ' e 


where is the constant pressure specific heat of the mixture of gases in 
the element. Upon insertion of Eq. (29) in Eq. (27) we find 


K 0 ~ K 0 

2 . 
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(30) 


Then, combining Eq. (26) and Eq. (28-30) there obtains, after changing 
units, the equation for temperature rise in the element, 


dT _ , T 1 ^ 2 1/2 

dt “ h .325 *f. *0 0 

P l 2. 


1 ^7- (T - T.) 
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(31) 


which is in terms of the initial concentrations in the wake, Kf^ and Kq^ 
The cuiistant has value 
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for which R is the specific gas constant of the mixture of gases in the 
element. Equation (31) has English system units. 

From Eq. (31), we find the time for the temperature in the element 
to reach T by integration, obtaining 


p .325 e 1/2 ^ 

™TTV*r 1 ' 

4 f. 0„ 

i 2. 


(33) 


where 
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In this expression, x = T/E, E -- 2. 2 x 10^ °R, and 
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(35) 


With T properly chosen, the value of t computed from Eq. (33) may be 
regarded as characteristic of the induction time which precedes very rapid 
reaction. For example, T can be chosen as that temperature for which 
the integrand in Eq* (34) has its minimum* This would correspond to the 


moment of maximum rate of temperature rise, (dT/dt) 


max’ 
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Some simplification of Eq. (34) is possible by noting ihat e^^ x / 
decreases with extreme rapidity. For typical values of Q , 8q , andx., • 
the integrand, therefore, becomes very small before Q^x - x.) or 
0O 2 ( X " x p become significant compared to unity. It suffices, in general, 
to approximate Eq. (34) with 





in which x can be chosen as virtually any value for which e^ x / x*^ 2 <K 
l/x* / l/2 

e 1 V x^ 7 . Equation (36) is then a universal function of x- only, whose 
values appear on Fig. 21. 

Each element of reactive gases in the near wake of the parent drop 

begins at a temperature T. (or x.) from Eq. (21), and has an induction time 

given by Eq. (33). It is clear from Fig. 21 that ^{x^ } and therefore t, is 

extremely sensitive to T., and this in turn is mainly a function of initial 

mixture ratio. The mixture ratio is not only spacially variable within the 

wake, but its distribution is also not duplicated In practice between parent 

drops that are subjected to identical conditions, as is apparent from photo- 

9 14 

graphs such as Fig. 22. In spite of this, measured ignition times ’ 

are reasonably reproducible, and so the sensitivity to the actual distribution 


of fuel vapor in the wake is apparently reduced by the mechanism of ignition 


that is involved. 



An ignition mechanism which is consistent with these observations is 
as follows. It is postulated that "explosive ignition" of the entire near wake 
can be traced to homogeneous reaction which, occurring in a local reactive 
element, releases sufficient energy to produce a shock wave capable of 
initiating detonation in the remainder of the wake. For the present purpose, 
a simplistic initiation requirement will be assigned; namely, the local 
pressure rise must be equivalent to a Mach 3 shock wave. 

The induction period of the homogeneous reaction was assumed to occur 
at constant pressure. However, once reaction rates become very large 
(i. e. T —T), the process is better approximated by a constant volume 
assumption. If most of the temperature rise occurs at these high reaction 
rates, the total temperature rise is then simply 


T 

max 



C 


V 


(37) 


which again assumes the fuel -lean condition, and C v is the constant volume 
specific heat of the reactive element mixture. It is readily shown that if 


p 

i= -j . — 1 * 

the minimum mass fraction of fuel vapor required to produce a specified 
value of £ is, from Eq. (37), 
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(38) 


K, 


T. C £ 
i v 


w< 


For the equivalent of a M g = 3 shock, £ ~ 9. 

It is easily shown that the initial value of C in the element is related to 
Kj, through 

C y = K f C [1 + rt(y f / , (39) 


where y ^ and y ■ are the ratios of specific heats of the fuel vapor and oxidizer, 
and a is defined by Eq. (22) . Combining Eq. (38) with Eq. (39) and (21) , 
and assuming yJy — 1 results in 

X X 

K f =1/(1+ <4 , (40) 

min 


where 





(41) 


Elements with k. < k. do not possess sufficient energy to produce 
1 main 

the pressure rise corresponding to £. while elements with k - > Kf will 

f man 

require a longer induction time, because of their reduced T-. By taking 
\ (k.j) into account for fuel rich mixtures, a corresponding maximum k^ 
can be found in similar fashion, but the fuel-lean element will always ignite 
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first, due to its higher Mittal temperature. Hence in Eq. (40) is 

min 

used to compute the minimum ehemicai induction time prior to wake 
detonation. 

It is noted that the mass fraction of fuel vapor corresponding to the 
stoichiometric condition is easily found to be 

A. <£ * 

*f* = r+ A. <j> * * 

r Q 


(42) 


in which $ o * is the stoichiometric fuel-oxygen mass ratio, and A. is the 

initial mass fraction of oxygen in the oxidizer gas. If x f > x * explo- 

x min f 

sive ignition is not possible; i. e. more energy is needed to initiate detona- 
tion than is available in any reactive element. Each element will proceed 
with local homogeneous reaction. This is also what should occur in ele- 
ments having x^ < Kf m ^ n j which pass through (dT/dt ^ prior to r, and 
explains the appearance of luminous turbulent regions in the wake, which 
occur before explosive ignition is observed, cf. Fig, 20. 


C. Results and Discussion 

The time to explosive wake ignition, measured from first contact by 
the incident shock with the parent drop, is prescribed to be 


t. 


t* + 


chm 


m 


in which t* is the time to reach peak frontal diameter; namely, 
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•*'*'■ •* > *r j' ■rTT'Z - l __ rL . . <t ^ __ , ^ ^ ^ ^ .-■,— 

t* = (D q T*)/(u 2 p l/2 ) , (44) 

a nd T c k m is obtained from Eq. (33) and (36) using Kf . from Eq. (42) 
to compute T\ in Eq. (21). The results from this computation are compared 
with experimental data on Fig. 23-25, showing the effects of Mach number, 
initial oxidizer pressure, and parent drop size. 

The phenomenological description given here is not applicable to shock 
wave interactions with parent drops under conditions which would produce 
first generation juspray whose characteristic quiescent evaporation timos 
are significantly less than t*. In such a case, explosive ignition 
within the recirculation zone, prior to t*, is conceptually possible. Under 
typical conditions, this roughly restricts the present formulation to D q > lOG/i. 
It should also be pointed out that the chemical induction time is quite sensitive 
to the liquid boiling temperature, T^, through its effect on T^, while accurate 
high temperature boiling point data is not available for all hydrocarbons. 
Moreover, the Edelman- Fortune reaction rate equation is strictly applicable 
to hydrocarbons (both olefin and parafin) , and to the approximate temperature 
range 1440 < T < 5400 °R. Other reaction rate expressions must be used 
to derive induction time expressions similar to Eq. (33) for other fuels and 
temperature ranges. 
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2. Self Oxidizing, Bi-Liquid (SOBL) Spray Detonations 


All of our earlier work on heterogeneous detonation was concerned 
with liquid fuel drops distributed in a gaseous oxidizer. In these studies 
it was found that such detonations are relatively easy to generate. The 
interest of this extension, then, was to investigate the possibility of 
detonation in a system comprised of liquid fuel drops, liquid oxidizer 
drops, and an inert gas. In such systems there exists the possibility 
of realizing much higher pressures, the limit being those associated 
with condensed detonations as the inert gas concentration approaches 
zero. In order that SOBL detonations be realized, it is required that 
the shock wave create a microspray in the wake of the fuel and oxidizer 
drops, that this spray evaporates and mixes with the other species, 
and, further, that all of this occurs rapidly enough that heat and momen- 
tum losses to the walls of the detonation tube don’t preclude ignition. 

The necessity of mixing the wakes of the fuel and oxidizer drops implies 
that the reaction zone length will be much longer, that wall effects will 
be more pronounced, and hence that detonation will be more difficult. 

Kerosene and concentrated hydrogen peroxide (50%) were chosen 
for the liquids with the gas being nitrogen. The vertical detonation tube, 
15/8 in. by 1 5/8 in. and 12 ft long, was used for these studies. Exten- 
sive modifications were made to the drop generator so that the desired 
drop sizes and amounts for the two different liquids could be accommodated. 


Further, a blast wave driver was added to the bottom of the vertical 
tube so that a strong blast wave could impact the heterogeneous mixture. 
The blast wave was generated by the detonation of and at initially 
elevated pressures. 

Hydrogen peroxide (50%) is a strong oxidizer which presented 
difficult storage and handling problems, particularly when used in the 
drop generating system. The dissociation of H^O^ sa ^ e ^ problems 

(and one accident) and also upset the requisite hydraulic character of the 
drop generating system. However, these problems were finally over- 
come and the experiments conducted. 

.An analysis was performed to guide the experiments and particular 
attention was given to the global equivalence ratio, the ratio of the 
diameter of fuel drops to that of the oxidizer drops, and the ratio of the 
maximum interaction distance of the two liquids to the hydraulic radius 
of the tube. While difficult to judge in some cases, the conclusion 
reached was that detonation was not realized. However, these studies 
were not all inclusive and it is believed that SOBL detonations can be 
generated under the appropriate conditions and further investigation is 
warranted. Stronger peroxide solutions, other liquid combinations, 
smaller drop sizes, and larger detonation tube cross section would be 
appropriate steps in this direction. A complete write up on this 
investigation is being prepared for publication and should be available 
in the near future. 
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APPENDIX A -—DROP GENERATOR POWER REQUIREMENTS 


1. THEORY 

A. Power and Amplitude of Disturbance 

Consider a jet of fluid emerging from a vertically mounted capillary 
tube in a coordinate system moving at the same velocity a s the jet. Prior 
to the application of any disturbance the fluid inside the stationary column 
(which is what the jet appears to be) has a pressure which exceeds the ambient 
pressure by an increment attributable to the surface tension. This incre- 
ment is inversely proportional to the jet radius. Hence if the pressure in 

the jet is made to vary in a regular manner by some means, variations in 

t 

jet radius and jet instability may result. Rayleigh showed that for a given 
jet diameter the disturbance wavelength for which the amplitude will grow 
most rapidly is given by, 

X R = 4. 508 d. ’ (I) 

where is the wavelength of the disturbance with maximum growth rate 

d. is the diameter of the jet 
The wavelength, frequency and jet velocity are related by, 

\ « u./f 

or from Eq. (1) 

j L = u./4, 508 <L 

it ) 3 
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( 2 ) 




where is the frequency for maximum disturbance growth rate, the 

Rayleigh frequency 
u. is the jet velocity 

Equating the volume of a drop to that contained in one wavelength 


l*n 3 _ £ \ a 2 
6 0 4 j 


The jet diameter may be replaced by the capillary diameter without appreci- 
able error to give 

D 0 /d t = 1. 89 (3) 


where Dq is the drop diameter 

dj. is the capillary (or needle) diameter 

For a jet of liquid of viscosity density p^, and surface tension cr injected 

into an atmosphere of viscosity fi and density p , an expression for th-3 

2 

minimum power required for jet breakup is, 


V . p * 

mU1 64 vl r t 2 I? 1 


(4) 


where (see Fig. A-l) the needle (capillary) is mounted in a fixed platform, 

Rg, and is connected to the supply tank through a supply tubing, part of 
which is mounted on a platform, of mass, m, being vibrated at frequency 1 
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l?min average minimum power required for jet breakup 

m mass of vibrated assembly 

r radius of vibrated tube 

L length of vibrated tube 

v ^ kinematic viscosity of the liquid 

p^ density of the liquid 

cr surface tension of the liquid 

f frequency of vibration 

q growth parameter in Rayleigh’s analysis 

x breakup distance from needle tip 

= 8x 
n 

x^ natural breakup distance 

e x/x 

n 

The natural breakup distance is the distance in which the jet would 
disintegrate naturally without a force disturbance being applied. Grant and 

3 

Middleman presented a comprehensive study and observed that for a fixed 

jet size d,, the breakup distance is a function of jet velocity u. and has a 
•* 1 

maximum at the point u. = u.*. For u. < u * 

F 3 3 3 3 

3 = 19. 5 (We 1 / 2 + 3 We/Re) °* 85 (5) 

t 

2 

where We is the Weber number = (p^ u d^/o) 

Re is the Reynolds number - (p^ u 


also m* corresponds to the condition, 

Re* - 325 (Z*) 


-0.28 


where Re* is the Reynolds number corresponding to u.* 

3 

l/2 

Z* is the Ohensorge number (We 7 /Re) corresponding to u.* 

3 

For u. > u.* or Re > Re*, x /d, cannot be evaluated in a closed form and 
3 3 n t 

2 

an iteration procedure has to be used . 

Since the power required depends on the mass of the assembly, which 
will depend on the kind of apparatus used, a more useful parameter is the 
amplitude of vibration, A, which is coupled to the power, P, through a 
simple relation for a sinusoidal disturbance, 


a 2 = 


q 2 .3 

87T mf 


( 6 ) 


combining Eqs. (4) and (6) 


a 2 -qx/u. 
(3<j r e 3 

^min ” 16 J 2 irft r^. Lf 


( 7 ) 


The pressure perturbations produced in the liquid as a result of the vibra- 
tions are given by 




16 V2 7T ^ ALf 


p ( 16 VI 7T pU ALf 

P 0 P Q P r 2 

where Pq is the atmospheric pressure. 
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B. Dimensional Analysis 

The important variables are d fc , u , and p ^ which will be used as 
repeated variables along with d ? L, A, o, andx. This leads to the 
dimensionless parameters 

,,T 1 = d/ d t y v 2~ » 73-g = A/d fc 

\ = ( u j/ ~ 4. 508 (for f - f R ) , 7Tg -• (P £ d fc u./ p^) - Re 

2 

TTg - (p^ d t U. /cr) a We , tt 7 a x/d fc 


Then from Eqs. (7) and (8), a straightforward substitution gives 

1/2 

/V*Y 

A-k ) 


( 10 ) 


and 


0 . 948 7T>- 


7T, 


3min 


K i W 

(^g) 


3/2 


-* e . 


u 6 ) 


172 


at) 
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where 


K 


if* 



1 64 ir 2 


(12) 


Now to study the dependence of amplitude on other variables, ideally 
it is necessary to vary one while holding all others constant. In the experi- 
ment reported here the ratios 7r^ and were kept constant (tt^ = 2, ^ = £57). 
Also, since it was desirable to operate at for any velocity u., tt ^ was also 
a constant at a value 4. 508. Then Eq. (11) predicts the dependence of 
7 Tg on 7r 5 , TTg, and 7iy • 


2. EXPERIMENTAL ARRANGEMENT 
A. Apparatus 

Apparatus reservoir R contained prefiltered test liquid. It was pressur- 
ized (usually to a constant 30 psig pressure) using a nitrogen bottle S. The 
supply manifold consisted of 1/ 4 in. ID clearflow tubing up to the coarse 
metering valve and 0.065 in. ID Teflon tubing from to the vibrated 
tube and from to the stationary capillary (needle) N^. 

The platform P^ was attached to a M. B. Manufacturing Co. (Model 
631H) shaker , V t which was excited by an audio oscillator (HP Model 2020) 
and an amplifier (MB, Model 1132034). The platform vibrations were 
sensed using a 6 mv/g sensitivity accelerometer, A (Mefcrexlnst. Co. 

Model 512-1) . The accelerometer output and the vibrator input current (£or 
testing its sinusoidality) were displayed on a dual beam oscilloscope, O. 


The pressure drop across was measured using a zero volume 
Micro Systems pressure transducer. The pressure transducer output was 
displayed on a Fluke 8300 A DVM. The sensitivity was 1.95 mv/psi. The 
breakup pattern was observed using a General Radio Co. Strobotac Type 
1531 -A. 

The shaker was mounted on isolation mounts to isolate the shaker and 
the rest of the system. The entire system was kept in a room maintained 
at 70°F to eliminate any temperature variation effects on viscosity and other 
parameters. 

B. Measurements and Calibration 

It was necessary to obtain a calibration of pressure drop across the 

needle (measured by transducer PT) against the average jet velocity. 

This was done by measuring the mass flux through the needle (for all test 

* 

liquids) for various pressure settings. Since the needle diameter was fixed 
and the measurements were to be made at the Rayleigh frequency*, this 
calibration could be converted to a P gefc vs frequency curve. 

Since the amplitudes of vibration involved were extremely small they 
were not measured directly. Instead the acceleration of the platform was 
measured using the accelerometer A. The output of the accelerometer 

*For a given frequency, various settings of pressure drop were tried 
and the one corresponding to minimum power was found to be approximately 
the same as predicted by Rayleigh. 
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was directly measured on the scope and the acceleration amplitude was 
converted to vibration amplitude knowing the frequency of vibration. They 
are coupled by 

a = 4 t? -f^ • A (13) 

where a is the amplitude of acceleration. 

The power measurements were made for various frequencies with the 
pressure drop (or velocity) corresponding to that which was found from the 
previously made operating curve. Measurements were made only up to the 
point beyond which either the power supplied by the amplifier or the vibra- 
tions (as depicted by the accelerometer output on the scope) were no longer 
sinusoidal. 

Once the flow rate was adjusted for a given frequency, the power was 
gradually increased until a uniform string of drops was produced. When 
the stroboscope is used as a source of light to view the breakup pattern, 
a uniformly broken jet produces a stationary pattern of drops as in Fig. A -2. 
When the power is slightly reduced the disturbance introduced is incapable 
of producing uniform drops, even though it still breaks up the jet as can be 
seen from Fig. A-3. When the power is completely shut off the jet breaks up 
"naturally" and since an ordinary lamp was used as a light source for this 
part, the long exposure time produces the apparent multi- exposure pattern 
shown in Fig. A -4. 
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The patterns such as in Fig. A-2, were taken as acceptable and the 
minimum power required to produce them was recorded. If the power is 
increased beyond this value no significant changes take place in the appear- 
ance of the pattern except that the point where the jet begins to breakup 
shifts towards the needle, i. e. the breakup distance x is reduced. 

It was also observed that the breakup distance, x, corresponding to 
the acceptable breakup pattern is always less than the natural breakup 
distance, x . However, these distances were not measured. 

3. CHOICE OF PARAMETERS AND RESULTS 
A. Choice of Liquids and Tube Sizes 


TABLE A-l. PHYSICAL PROPERTIES OF LIQUIDS TESTED* 


No. 

Substance 

Density 

Slugs/ft 3 

Viscosity 

lb-sec/ft^ 

Surface Tensi< 
lb/ft 

1 

Kerosene 

1.44 

3.44 xl(T 5 

1.94x 10" 3 

2 

Ethyl Alcohol 

1. 52 

2. 50 x 10"'* 

1.50xl0“ 3 

3 

Benzene 

1. 74 

1.31 x 10" 5 

1.98 x 10" 3 

4 

Water 

1.94 

2.09 x 10" 5 

5.01 x!0“ 3 


^Actual measurements for the sample were used for kerosene. Otherwise 
the source was the Handbook of Chemistry and Physics. 

The liquids were chosen so as to have pairs in which the two liquids had all 
properties approximately the same except one. Thus, for viscosity, the pair is 
benzene and Kerosene, and for surface tension, alcohol and water. 
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The power is extremely sensitive to the needle diameter (see Ref. 2, 
p. 236) and hence a needle diameter of 0.027 in., which produced 1300 jmm 
drops (see Eq. (3)) and which required a reasonable amount of power for the 
chosen frequency range, was selected. 

For all liquids except water the tests were carried out at this needle 
diameter. For water an additional needle of 0.039 in. diameter was used. 
The length of the needles was 2 in and they were made by telescoping two 
or three stainless steel tube pieces to give an OD of 0. 016 in. to fit the 
supply tubing and the required ID. The vibrated tube was mounted on 
(Fig. A-l) and had an ID of 0. 054 in. and was 1 in. in length. The supply 
tubing downstream of was flexible enough so as not to resist the vibra- 
tions. The entire fuel system was maintained hydraulic so as to obtain 
instantaneous response to any changes in flowrate (using either C % or C^) . 

B. Results 

From the recorded data, values of the various parameters, ?r^ - 7Tg, 
were calculated. Because of the choice of needle diameter and Rayleigh 
frequency, and i were constant throughout except for water. 

Thus, 

7T^ = 2. 0 , 7r 2 = 37. 0 , = 4. 508 

for all liquids. In addition, for water, 

7r x = 1.385 , ^ - 25. 6 , 7^ = 4.508 
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From measured values of parameters - tTq, the values of were calcu 
lated using the relation: 



<V 


1/2 


0.948 


log 


K 1 ^ 


3/2 


V*3 


(14) 


which is just an inversion of Eq. (11). 

Values of dp were calculated using Eq. (5), or iteration procedure, 
as appropriate for the same conditions knowing the Re and We (i. e. 7 Tg 
and Tig) and from this values of 9 = (x/x^) were obtained. 

Also, the values of P’/Pq were calculated from Eq. (9), using from 
the measured or known values. These values of tt^, tTq, (x /dp, 

0, and P t /Pq are tabulated in Table A-2. 

Comparing the values of and (x / dp it can be easily seen that 
n'rj < ^/^p. This is an expected result, since the externally introduced 
disturbance should break up the jet earlier than if it were broken up naturally. 
In fact, this is the reason why any external power is needed at all. This is 
in agreement with the previously recorded observation that increase of 
power after the minimum point merely results in reducing the breakup 
distance. The regularity of the disturbance results in a regular breakup 
pattern. 

Thus in Eq. (4), x can be replaced by 0 • x^, which is a modification of 
the original equation from Kef. 2. Note that the only way to substitute a 
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meaningful distance x from a purely theoretical basis is by relating it to 
the natural breakup distance. Thus, the original equation in Ref. 2 could 
be looked upon as a zeroth order approximately where the a priori assump- 
tion, x = x^, was made. However, as explained above, it would be expected 
that (x/x n ) < 1, always. 

It can also be seen from Table A-2 that the value of 9 varies between 
0.4 and 0. 8, being very close to 0. 8 for many readings. Therefore, the 
average value of 9 - 0. 58 was chosen to relate x with x . This could be 
taken as a first approximation. With this x the values of ir^ ~ A/d^ were 
calculated theoretically and are tabulated in the last column of Table A-2. 
The experimental values are compared with theory in Figs. 5 and 6 for 
the four liquids tested. The disagreement for some points is considerable, 
particularly when 9 takes extreme values like 0.4 and 0. 8. This is due to 
the exponential dependence of on n ^ (i. e. 6 ) . From the values of P t /Pq, 
appearing in Table A-2, it can be seen that the pressure perturbations 
produced are indeed very small and are comparable with those in an audible 
sound wave. This justifies the small perturbation assumptions made in 
deriving Eq. (9), and others. 

4. DISCUSSION 

It seems that a better correlation between the actual breakup distance x 
and the natural one x^, i. e. the value of parameter 0 and its dependence on 
flow parameters like Re and We determined from actual experimental 
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measurements, would improve the accuracy of the theoretical predictions 
of required minimum amplitude (or power) considerably. It is believed 
that this uncertainty is the largest contributor to the discrepancy between 
experiment and theory. Further investigation is required. 

It should also be noted that the analysis presented here could be used 
irrespective of the method of introducing the disturbances, provided they 
are introduced axially. 


TABLE A-2. TABULATED RESULTS 


^3 

ff 5 

*6 

*7 

x/ d f 

n l 

9 = x/x 

n 

p, / p o 

Theo. 

^3 

x 10 " 4 
0.638 

576.0 

61.9 

64.94 

130.6 

0. 50 

x 10' 6 
0. 323 

x 10' 4 
0.168 

0.749 

672.0 

84.3 

73.48 

141.9 

0. 52 

0.409 

0.292 

0.720 

768.0 

110.1 

83.67 

149.4 

0. 56 

0.421 

0. 538 

0.797 

864.0 

143.2 

93.01 

153.0 

0.61 

0. 507 

1.095 

0.645 

960.0 

172.0 

104. 56 

153.1 

0.68 

0.421 

1.972 

1.30 

1056.0 

208.0 

103.66 

152. 5 

0.68 

0.886 

3.418 

1.79 

1152.0 

247.6 

107.05 

151.9 

0.70 

1.276 

5.49 

2.31 

1248.0 

290.7 

110. 59 

131.5 

0.71 

1.724 

8.04 

4.23 

1342. 5 

337.1 

106.68 

151.1 

0.70 

3.272 

10.97 

4.70 

1440.0 

387.0 

111.42 

150.7 

0.73 

3.758 

14.74 

6.48 

1536.0 

440.3 

111.04 

150.4 

0.73 

5.348 

18.67 

10.4 

920.0 

37.3 

40.22 

108.4 

0.38 

1.222 

0.303 

4.87 

1380.0 

83.9 

65. 70 

134.8 

0.49 

0.700 

1.31 

6. 57 

1670.0 

114.2 

72. 39 

138.9 

0. 52 

1.021 

3.12 

7, 61 

1840.0 

149.1 

79.97 

138.6 

0. 58 

1.264 

7.25 

6.48 

2070.0 

118.1 

77. 72 

138.0 

0. 56 

0.713 

5.17 

6.40 

2300.0 

233.0 

100.97 

137.6 

0.73 

1.187 

23.4 

6.69 

2530.0 

281.9 

109.43 

137.2 

0.79 

1.302 

35.6 

7.99 

2760.0 

335. 5 

115.12 

136.9 

0.80 

1.623 

49.9 

10.89 

844.0 

14.4 

28.10 

74.0 

0.38 

1.861 

0.263 

17.46 

1055.0 

22. 5 

32.20 

85.0 

0.38 

3.337 

0. 561 

17.85 

1266.0 

32.4 

37.96 

92.0 

0.41 

3.737 

1.34 

10.79 

1688.0 

57.6 

53.49 

94.9 

0. 56 

2.608 

8.61 

16.97 

1879.0 

72.9 

55.43 

94.7 

0. 58 

4.420 

17.6 

168.7 

3165.0 

202. 5 

54.30 

94.2 

0.58 

55. 837 

165.4 

416.0 

4220.0 

360.0 

51. 46 

94.0 

0.55 

159.0 

357.4 

0.40 

906.0 

18.3 

45.98 

79.0 

0.58 

0.031 

0.416 

44.7 

1208.0 

32. 5 

32.05 

90.1 

0.36 

4.01 

2.12 

25.7 

15 \0 

50.8 

43.40 

92.3 

C. 47 

2. 578 

6. 53 

15.2 

1812.0 

73.1 

55.98 

92.7 

0.60 

1.669 

19.14 
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Source 


Figure 1. Needle-Electrode Assembly. 
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Figure 2. Needle Tip Configuration. 
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Figure 3 e Equipment Arrangement. 
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Figure 4. Shock Wave Interaction of Initially Ignited Drops, 
P = 30 in. Hg 0 0 , M = 2.02, d =1.6 mm. 
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Figure 5. Streak Schlieren Photograph of Non-ignited Drop 

V, = 30 in. Hg 0 0 , M o = 1.73. 


iression Wave 


Figure 6. Streak Schlieren Photograph of Drop with Wake Flame 

P 1 = 30 in. Hg 0 9 , M c = 1.72. 


liu ul( 111 

Reference 


Mark on Window 














d/d 



0 20 40 60 80 100 


At C/i SI 

Figure 7. Drop Deformation, 
Wake Flame Conditions. 



Figure 8. Drop Displacement Parallel 
to the Direction of Flow. 
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Figure 9. Limiting Initial Pressure of Pure Oxygen for Blow Off. 
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Figure 10. Comparison of Initially Non-Ignited and Ignited Drops at 
M = 2.05, P = 29. 1 in. Hg in Nitrogen, Air and Oxygen. 
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Figure 10. Concluded. 
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Figure 11. Strong Shock Wave over Non-Ignited Drops in N ? and 0 ? ; 

M = 3. 5, P i = 29. 1 in. Kg. 
s ’ 1 b 


72 







' Hnf ' ' 



E|| 

wm 





(a) 


(e) 


Figure 11. Concluded. 
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Wake Shock 


Figure 13. Shock Wave over a 

Water Drop, d = 750 (im, 

M = 2. 7, t = 2. 6 us. 
s 


Figure 14. Shock Wave over a 
Water Drop, d = 750 pm, 

M g = 2. 0, t = 15. 8 ps. 
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(b) 

Figure 15. Shock Wave over a Water Drop, 

(a) d = 750 pm, M s = 2. 7 , t = 4.4 ps 

(b) d = 1400 pm, M s - 3. 34. 
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Figure 20. Detonation of Diethylcyclohexane Drops i„ Oxygen, d = 


2600 /urn. 
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Figure 22. Shock Wave over Water Drops 
M = 3. 25, t = 38. 8 jus. 















Figure A-l. Schematic of the Apparatus 
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Figure A -2. Drop For nation with Figure A -3. Drop Formation with Figure A-4. Natural Drop 

Correct Power Setting. Insufficient Power. Formation — No Power. 
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Figure A-6. Comparison of Experimental and Theoretical Results 

Ethyl Alcohol, Benzene. 


89 



